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Abstract: A recent study bv lWvlie et al.l (|2006T l has revealed that s-process element abundances are 
enhanced relative to iron in both red giant branch and asymptotic giant branch stars of 47 Tucanae. 
A more detailed investigation into s-process element abundances throughout the colour-magnitude 
diagram of 47 Tucanae is vital in order to determine whether the observed enhancements are intrinsic 
to the cluster. This paper explores this possibility through observational and theoretical means. The 
visibility of s- and r-process element lines in synthetic spectra of giant and dwarf stars throughout 
the colour magnitude diagram of 47 Tucansa has been explored. It was determined that a resolving 
power of 10 000 was sufficient to observe s-process element abundance variations in globular cluster 
giant branch stars. These synthetic results were compared with the spectra of eleven 47 Tucanae giant 
branch stars observed during the performance verification of the Robert Stobie Spectrograph on the 
Southern African Large Telescope. Three s-process elements, Zr, Ba, Nd, and one r-process element, 
Eu, were investigated. No abundance variations were found such that [X/Fe] = 0.0 ± 0.5 dex. It was 
concluded that this resolving power, R ~ 5 000, was not sufficient to obtain exact abundances but 
upper limits on the s-process element abundances could be determined. 

Keywords: instrumentation: spectrographs — techniques: spectroscopic — stars: abundances — 
globular clusters: individual (47 Tuc) 

1 Based on observations made with the Southern African Large Telescope (SALT) 
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1 Introduction 

Globular clusters are a rich area of study. They are 
very old and therefore can be studied as remnants of 
the early Galaxy and Universe. Their galactic orbits as 
individual clusters, and as a system of clusters, provide 
information regarding galaxy formation and structure. 
In addition, based on the assumption that all the stars 
in a globular cluster (GC) were formed from the same 
gas cloud, they provide a unique view of stellar evolu- 
tion and nucleosynthesis. Variations in chemical abun- 
dances between the stars can be investigated as evi- 
dence of the different types of nuclear processing that 
occur as the stars evolve, or as evidence of inhomo- 
geneities in the initial gas cloud. It is this latter area 
that the work presented here addresses. 

Extensive literature is available on abundance vari- 
ations, particularly c arbon, nitrogen and o xygen abun- 
dances, within GCs (jGratton et al.ll2004h . As larger 
telescopes and more sensitive instruments become avail- 
able, heavier elements have been added to the list of 
observed anomalies. Most recently, s-process element 
enhancements have been observed in asymptotic gi- 
ant branch (AGB) and red giant branch ( RGB) stars 
in 47 Tucanae (47 Tuc) (jWvlie et all 120061 ). For AGB 
stars with mass greater than 1 Mq the third dredge- 



up (TDU) occurs with the mixing of s-process elements 
formed in nuclear reactio ns in the C 13 pock et into the 
observable surface layers (jBusso et al]|200'lT l . However 
the AGB stars in 47 Tuc have mass less than 1 M©, 
making this an unlikely caus e of the observed enhance- 
ments (jGilliland et al] 1 1998ft . There is no equivalent 
mechanism that would produce the observed enhance- 
ments in the RGB stars. The most likely explanation 
for the enhancements is that the gas cloud from which 
the stars formed was itself s-process enhanced. Al- 
ternatively, t he stellar atmospher es were polluted by 
stellar winds (|Cannon et al.lfl998T ). 

Until relatively recently high resolution spectra could 
only be obtained for the brightest giants in GCs as the 
remaining stars were too faint for this type of obser- 
vation. As the bright giants are evolved stars their 
atmospheres may have had the products of internal 
nuclear burning mixed through them, modifying their 
spectra. However, stars on the main sequence (MS) 
have not yet experienced mixing and it can be as- 
sumed that their spectra reflect the chemical compo- 
sition of the initial gas cloud. 10 m class telescopes 
have the capabilities to observe these fainter dwarf 
stars in GCs. By comparing the spectra of main se- 
quence stars to giant branch stars in the same cluster 
the intrinsic chemical composition of the cluster can 
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be investigated. A new technique which pursues this 
comparison is to analyse the integrated light from a 
GC. This is a technique that is being developed and 
will a lso allow comparison of galacti c to extra-galactic 
GCs (|Mc William fc Bernstein]|2007ft . 

The purpose of the current study has been to mo- 
tivate a medium resolution study of 47 Tuc using the 
multi-object spectroscopy mode of the Robert Stobie 
Spectrograph on the Southern African Large Telescope 
and the AAOmega spectrograph on the Anglo- Australian 
Telescope. 

2 Telescope & Instrument Combinations 

47 Tuc is the second largest observable galactic globu- 
lar cluster and is only visible in the southern sky. Both 
the Southern African Large Telescope (SALT) and the 
Anglo- Australian Telescope (AAT) are ideally situated 
to perform a comprehensive survey of s-process abun- 
dances in this cluster. 

SALT is a 10 m class telescope with instru menta- 
tion a t the forefront of astronomical technology (|Bucklev et 
l2006ft . Currently in operation on SALT is the Robert 
Stobie Spectrograph ( RSS) which has low t o medium 
resolution capabilities (|Nordsieck et aLl200lh . The high 
resolution spectrograph (HRS) has completed the de- 
sign stage (|Cottrell et al.l [20051 ) and construction has 
commenced at the University of Durham in the UK. 

AAT, a 3.9 m class telescope, has capabiliti es for 
mult iple object surveys using AAOmega and 2dF (ISharp et 
I2006T ). These instruments have already been used suc- 
cessfully for element abu ndance surveys of g lobular 
clusters, including 47 Tuc (jCannon et al.ll2003f ). 

This paper considers the feasibility of using RSS 
and AAOmega for the 47 Tuc s-process abundance 
survey. These instruments have comparable maximum 
resolving powers, RSS at R ~ 10 000 and AAOmega 
at R ~ 8 000. Comparison will also be made with 
the high resolving power capabilities of SALT HRS, 
which has a nominal maximum resolving power of R ~ 
60 000. Section [3] addresses the question of whether or 
not variations in s-process abundances of at least +0.5 
dex will be visible in spectra obtained at the maxi- 
mum operating resolving power of these instruments. 
The line strength visibility for a range of key r- and s- 
process species will be explored in the effective temper- 
ature (T e ff) - surface gravity (log g) space represented 
by the colour magnitude diagram (CMD) of 47 Tuc 
for both high and medium resolving powers. Section|4] 
presents the results of abundance analyses of s-process 
elements in eleven giant branch stars in 47 Tuc that 
were observed using RSS on SALT in 2006. 

3 Stellar Atmospheres & Line Strength 

The chemical composition of a star is determined by 
the analysis of absorption lines in the stellar spectrum. 
Every line corresponds to the absorption of energy by 
an atom and there is a direct relation between the 
strength of a line and how much of the correspond- 
ing element is present in the star's atmosphere. How- 
ever there are other factors which affect the observed 
strength of a line. These include the temperature, 



pressure and opacity in the star's atmosphere. When 
considering the stars in a globular cluster the two key 
parameters for classifying the evolutionary stage of 
each star are the effective temperature, T e ff, and the 
gravity, log g. 



3.1 



Line Strength Analysis in the 47 Tuc 

'T e ff — log cf Space 



Figure [J is the CMD of 47 Tuc (| Cannon et al.lfl998l : 
lLelll977l ; Irlartwick fc Hesserjll974T )~ The stars range 
in T eff from - 4000 K to ~ 6000 K, and in log g from 
~ 1.0 to ~ 4.5. 



Figure 1: CMD of 47 Tuc l|Cannon et al.lll998l : lLeelll977t 

Hartwick & Hcsscr 1971) showing T e jj and log g (o) at 
which synthesised spectra are computed and stars observed 
using SALT RSS during performance verification (*). 
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The six circular points on the CMD shown in Fig- 
ure [1] were selected to represent key stages of stellar 
evolution at which s-process element abundance anal- 
ysis would be undertaken. The T e g and log g values 
for each point are displayed. The eleven star-shaped 
points are the locations of the eleven 47 Tuc stars ob- 
served during the performance verification of RSS on 
SALT which will be discussed in Section fl] 

The models and linelis ts used in the spectrum syn- 
thesis programme (MOOG lSnedenl (jl973l)) were sourced 
from the Kurucz website ( |http: / /kurucz. harvard. ediij ) . 
Each model had a metallicity of [Fe/H] = -0.5 which 
is close to th e metallicity of 47 Tuc ([Fe/H] = -0.7) 
((Harris! [l996f l. A microturbulence of 2.0 kms -1 was 
selected. The Kurucz linelists are semi-empirical and 
some adjustment of oscillator strength (log gf) was 
requi red for individual lines (|Kurucz fc Pevtremannl 
Il975f ). In particular, they were adjusted to include 



ements dBiemont et al.l 


19811; iDen Hartog et al.ll2003l; 


lHannaford et all Il982l: 


Lawler et all 1200 if). Initially 



the linelists and models were used to produce spec- 
tra that were compared to the high resolution atlas 
of Arcturus which is available on the NOAO website 
(http://www.archive.noao.edu). Arcturus has a simi- 
lar metallicity to 47 Tuc stars. The log gf values of 
other species in the linelists were adjusted to provide 
a best fit to Arcturus. 
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Synthetic spectra were computed for each identi- 
fied evolutionary point in Figure [T] Three spectral re- 
Ions were synthesised based on lines used in lWvlie et al.l 
20061 1. The key species present in each of the spectral 
regions are set out in Table [T] where light s- represents 
the elements Y, Sr and Zr, and heavy s- represents Nd, 
La and Ba. 



Table 1: 


Key features in 


synthesised 


regions. 


Region (A) 


Feature (A) 


Species 


Process 


6139 - 6145 


6140.46 


Zr I 


light s- 




6143.18 


Zr I 


light s- 


5803 - 5807 


5804.00 


Nd II 


heavy s- 




5805.77 


La II 


heavy s- 


6644 - 6646 


6645.13 


Eu II 


r- 



Model spectra were synthesised at high (R ~ 60 000) 
and medium (R ~ 10 000) resolving powers. The high 
resolving power provides comparis on with a current 
instrument (UCLES) at the AAT (jWvlie et al.l l2006h 
and an equivalent instrument, SALT H RS, currently 
under contruction (jBucklev et al.l[2006l ). The medium 
resolving power is about the highest achievable with 
RSS on SALT and comparable to AAOmega's maxi- 
mum resolving power on the AAT. 

Each of Figures [5] to fl] and Figures [5] to [TT] corre- 
sponds to one of the six stellar models shown in the 
CMD of 47 Tuc in Figure [1] In each spectral region the 
abundance of the key species was varied with respect 
to the model abundance. 

By comparing the spectral regions between the mod- 
els, and between the two resolving powers, the changes 
in line strength for the different species are very clear. 

3.2 T e ff & log g effects on the GB 

Figures [2] [3] and |4] show the synthesised spectra for the 
three giant branch evolutionary stages. The progres- 
sion from the tip down the giant branch involves an 
increase in both T c g and log g. 

Both the neutral and ionised lines show a decrease 
in line strength down the giant branch. However the 
ionised lines are much less affected than the neutral 
lines. This is expected as the ionised state becomes 
dominant due to the increasing temperature. However 
the increased temperature also increases the H~ con- 
tinuous opacity which weakens both line species. 

The increased pressure has more complicated ef- 
fects as the neutral dominated populations are chang- 
ing to ionised dominated populations due to the in- 
creased temperature. Initially the neutral lines will 
weaken with the increase in pressure, but as the pop- 
ulation becomes ionised the effect becomes negligible. 
The ionised lines will initially weaken considerably with 
the increased pressure, but as the population becomes 
ionised this effect will lessen. 

The key observable in this sequence is the weak- 
ening of Zr I lines as the temperature and pressure 
increase, while the La II, Nd II and Eu II lines are less 
affected. The medium resolving power synthesis repli- 
cates the changes in line strength outlined for the high 
resolving power. All the features are identifiable and 



Figure 2: Synthetic giant branch spectra for key features 
computed at T e g = 4000 K, log g = 1.0. The dash-dot line 
represents a variation of —8.0 dcx simulating when none of 
the species, X, is present. The solid line represents no varia- 
tion of the species abundance from the model ([X/Fe]=0.0). 
The dashed line represents a variation of +0.5 dex and the 
dotted line, +1.0 dex. 
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Figure 3: As for Figurc[2]but T eff = 4500 K, log g = 1.5. 
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Figure 4: As for Figure[2]but T off = 5000 K, log g = 2.0 
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distinguishable at the +0.5 dex level, although the Zr I 
lines at the highest temperature have weakened con- 
siderably and would not be of use as Zr abundance 
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indicators at medium resolving power. 



4 SALT Performance Verification 

In June and October of 2006 spectra for eleven giant 
branch stars in 47 Tuc were obtained during the per- 
formance verification commissioning stage of RSS on 
SALT. Figure [T] shows the location of the eleven stars 
on the 47 Tuc CMD. The observations were carried out 
at two different camera settings obtaining two overlap- 
ping wavelength regions: 6020-6860 A & 5200-6150 A. 

These regions contain the key s- and r-process lines 
that were used in the study of theoretical line strength. 
The resolving power obtained for these spectra was ~ 
5000, which is not the maximum that will be possible 
with RSS. Signal to noise ratios between 50 and 140 
were obtained. 

Table [5] presents the effective temperatures and 
surface gravities for each star that were calculated from 
the stellar apparent magnitudes and colour ind i ces us - 
ing the calibration equations from lAlonso et al. 



Tabic 2: T e g and log g values from lAlonso et al.l (|1999T ) 
for each of the SALT PV Stars. 



Lee No. 



V 



B-V T off (K) 



3512 


11.79 


1.63 


4080 


0.7 


1513 


12.41 


1.32 


4450 


1.3 


2525 


12.43 


1.29 


4490 


1.4 


6519 


12.81 


1.27 


4520 


1.5 


6524 


13.08 


1.21 


4610 


1.7 


1506 


13.27 


1.15 


4700 


1.9 


3510 


13.63 


1.09 


4800 


2.1 


2604 


13.07 


1.04 


4890 


1.9 


4514 


14.35 


0.98 


5000 


2.5 


4515 


14.49 


0.98 


5000 


2.5 


4513 


14.59 


0.98 


5000 


2.6 



4.1 Temperature Sensitivity 

Figure [5] shows both the observed spectra of the eleven 
stars and the synthesised spectra of the region contain- 
ing key Zr I and Ba II features. 

The spectra are presented with the temperature in- 
creasing down the panel, replicating the traverse down 
the giant branch. This region contains two Zr I fea- 
tures and one Ba II feature. The vertical lines trace 
the line locations in each spectrum. 

In the coolest stars the neutral Zr features are quite 
distinct and get weaker until a temperature of around 
5000 K. The ionised Ba, while it does weaken slightly, 
remains apparent through to the hotter temperatures. 

This nicely replicates the theoretical analysis, show- 
ing the more dramatic weakening of neutral lines com- 
pared with ionised lines. It also is a good illustration 
of temperature sensitivity, the strength of the Zr I fea- 
tures providing a tem perature scale that c onfirms the 
scale derived from the lAlonso et al] (|l999f ) equations. 
The temperature sensitivity is replicated in the syn- 
thetic spectra. 



Figure 5: Temperature sensitivity: Zr I & Ba II. Key 
regions for observed (left panel) & synthesised (right panel) 
spectra. The key for the synthesised spectra as per Figurc|2] 
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4.2 Abundance Results 

Figures [6] and [7] present an abundance analysis of the 
two key spectral regions. In Figure the synthetic 
spectra (as for Figure[2j are overlaid onto the observed 
spectra for the Zr I & Ba II region for four of the stars. 
The error bars give an indication of the uncertainty due 
to noise in these spectra. 

Figure 6: Abundances: Zr I & Ba II. As for Figure[5]but 
observed spectrum overlaid on synthesised spectra. 
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Wavelength (A) 

The coolest star provides the best estimate. The 
Zr I features are prominent, again traced out by the 
vertical line. This spectrum best fits the model synthe- 
sis, indicating no Zr enhancement. The Ba II feature 
was also fitted best by the model abundance synthesis. 
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For the next two spectra, both of ~ 4500 K, while 
within error [Ba/Fe], [Zr/Fe] = 0.0 there is the possibil- 
ity of some Zr enhancement at a +0.5 dex level. 

For the hotter stars the Zr I features have dimin- 
ished too much to observe abundance variations. The 
Ba II feature is still prominent but the noise levels 
of these spectra do not give a convincing argument 
for anything other than model abundance. Overall an 
[Zr/Fc] = 0.0 with a possible uncertainty to + 0.5dex, 
while [Ba/Fc] = 0.0 for all stars. 

Figure [7] presents a similar abundance analysis of 
the region containing key Nd II & Eu II features. Even 
for the coolest star the Nd II features are too weak to 
provide a reliable analysis. The Eu II feature is more 
promising for the cooler stars (4500 K) with a good fit 
to the model abundance ([Eu/Fe] =0.0) within error. 
As the stars get hotter this feature also diminishes. 

Figure 7: Abundances: Nd II & Eu II. Key as for Figure[5] 

but observed spectrum overlaid on synthesised spectra. The 
observed spectra with T c g = 4890 K shows an absorption 
artifact due to a cosmic ray. 
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The analysis of these spectra shows that [X/Fe] = 0.0 
is the best solution in all cases with upper limits on 
the abundance values of +0.5 dex. For future RSS ob- 
servations obtaining the expected maximum resolving 
power of 10 000 will provide the required resolution to 
refine this upper limit. 

5 Future work - Main Sequence 

Comparison between the AGB and the RGB of a GC 
will provide a test for TDU if the AGB stars are s- 
process enhanced while the RGB stars are not. How- 
ever, as the spectra of GC MS stars correspond to the 
chemical composition of the initial GC gas cloud, to 
really investigate pollution versus nucleosynthesis the- 
ories the investigation must inevitably be extended to 
include the MS. The 10 m class of telescopes have the 
capability of observing these faint dwarf stars in dis- 
tant GCs. As discussed earlier, this investigation is 
also being pursued via integrated light techniques o n 
2.5 m class telescopes (|McWilliam fc Bernsteinll2007l '). 



Figure [8] shows the apparent magnitude limits of 
current (solid line) and future (dashed line) samples of 
47 Tuc stars from high and medium resolution instru- 
ments on SALT and the A AT. The SALT RSS limiting 
magnitude is expected to be ~ 22 mag which is suffi- 
cient for the spectra of main sequence stars in 47 Tuc 
to be obtained. 

Figure 8: V magnitude limits on 47 Tuc CMD with cur- 
rent (solid lines) and future (dashed lines) observations on 
SALT and the A AT. 
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5.1 T eff & log g effects on the MS 

Of the six evolutionary points on the 47 Tuc CMD in 
Figure [U one lies on the main sequence turn off, and 
two lie on the MS. Synthetic spectra were generated 
for these points as for the GB points in order to de- 
termine the visibility of s- and r-process lines at these 
temperatures and gravities. These spectra are shown 
in Figures I§1 to ITT1 

Progressing from the giant branch, the increased 
temperature causes increased ionisation which strength- 
ens the ionised lines. However the lines are weakened 
due to the increase in continuous opacity and the in- 
creased pressure to which ionised lines are sensitive. 
Figure [9] shows that the line strength of the ionised 
lines has reduced significantly in the high and medium 
resolving power spectra. 

The effect on the neutral lines is much more dra- 
matic. Both the increased ionisation and increased 
continuous opacity weaken the neutral lines but there 
is no effect at work to strengthen the neutral lines. Fig- 
ure [5] shows that in the high resolving power spectra 
the neutral lines have almost disappeared while in the 
medium resolving power spectra they are non-existent. 

Figures [TU] and [TJJ show the synthesised spectra 
for the two points on the MS. The models have the 
same log g value, so there is no change in pressure, but 
the T e g is lower for the point further down the MS. 
The differences in line strength between the models 
is due to the decrease in temperature. However the 
pressure is very high and the change from an ionised 
to neutral dominated population increases the pressure 
sensitivity of both the ionised lines and neutral lines. 

The neutral lines start to experience the weakening 
effects of the high pressure as the population returns to 
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being neutral dominated. However the decreased tem- 
perature means decreased ionisation and a decrease in 
continuous opacity. 

Figure 9: Synthetic Spectra for key features computed at 
T e jj = 5500 K, log g = 3.5 for a main sequence turn off star 
with key as for Figure [2] 
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Fi gure 10 '. As for Figure [9] except on the main sequence 
at T eff = 5000 K, log g = 4.5. 
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Figure 11: As for Figure [9] except on main sequence at 
T eff = 4500 K, log g = 4.5. 
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These act to significantly strengthen the neutral 
lines on the descent down the MS (see Figures [10] and 
llll l. The Zr I features have become very pronounced 



in the lowest temperature model. For the ionised lines 
the comparison between Figures QUI and fTTI shows that 
overall there is some strengthening of the ionised lines 
due to the decrease in temperature but the strong pres- 
sure ensures it is a small effect. 

This comparison indicates that neutral lines are 
the best candidates for abundance analysis of MS stars 
in both the high and medium resolving powers. The 
resolving power of SALT HRS (R ~ 60 000) will be 
required to fully test the hypothesis of intrinsic MS en- 
hancements using these features. However those MS 
stars above RSS's limiting magnitude but still suffi- 
ciently low in temperature will have neutral lines that 
are strong enough to be observed at medium resolving 
power. 



6 Conclusions 

The analysis of these initial (PV) SALT RSS observa- 
tions of eleven stars in 47 Tuc indicate no enhancement 
of Zr, Ba or Eu, that is [X/Fe] = 0.0 for these three 
species. However, the uncertainty is sufficiently large 
(+0.5 dex) that there may be some star-to-star varia- 
tion. This will have to await the more extensive survey 
at the maximum resolving power (~ 10 000) of RSS. 
The SALT RSS PV observations do provide an excel- 
lent temperature sequence and co nfirmation of the ef - 
fective temperatures derived from lAlonso et al.l (|l999h 
calibrations. 

The exploration of the T e g — log g space in terms of 
line strength has significantly refined the line selection 
process for the abundance survey. The analysis shows 
that neutral lines are observable over a greater range 
of log g values than the ionised lines, but within a 
smaller range of T e g values. Conversely the ionised 
lines are observable over a greater range of T e f{ values 
than the neutral lines, but within a smaller range of 
log g values. This broadens the range of temperatures 
and gravities from which the survey star sample will 
be selected. 

The characteristics of the instruments under con- 
sideration further limit the star sample. The medium 
resolving power of RSS and the magnitude limit of 
SALT do reduce the sample options with regards to 
the MS. There is a small range on the MS where the 
temperatures are cool enough and the magnitudes are 
bright enough that abundances could be obtained us- 
ing RSS, though SALT HRS is more suited to ob- 
taining this data. However, s-process element abun- 
dance analysis to just below the horizontal branch, 
which this analysis shows is quite possible using RSS or 
AAOmega, will provide definitive results on the nature 
of s-process element abundance variations in 47 Tuc. 

The next stage of this investigation is to consider 
globular clusters over a range of metallicity. The stel- 
lar models considered so far have a metallicity similar 
to that of 47 Tuc. The s-process element abundance 
survey results for 47 Tuc will then need to be compared 
to observational surveys of other clusters at different 
metallicity. This will show whether any s-process ele- 
ment abundance variations are isolated to 47 Tuc, or 
are a general characteristic of globular clusters. 
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